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A zone of reactions is determined and then exploited as a tool in studying the space-time structure 
of an interacting system formed in a collision of relativistic nuclei. The time dependence of the 
reaction rates integrated over spatial coordinates is also considered. Evaluations are made with the 
help of the microscopic transport model UrQMD. The relation of the boundaries of different zones 
of reactions and the hypersurfaces of sharp chemical and kinetic freeze-outs is discussed. 
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Introduction. In the collision of nuclei at high en- 
ergies, a strongly excited system of interacting particles 
is formed. In fact, a fireball is identified with the zone 
of reactions by such a definition, i.e. with a space-time 
region, in which the reactions of particles occur. Hence, 
the zone of reactions must reflect the space-time char- 
acteristics of a fireball, and its study gives information 
about the evolution of the interacting system. 

While studying the evolution of a fireball, it is impor- 
tant to know the size of the regions where the majority 
of the various processes are running. Depending on the 
model describing the system, we can distinguish the re- 
gions of the formation of a fireball, its isotropization and 
thermalization, the creation of particles, the regions of 
a chemical freeze-out and a kinetic one, etc. This al- 
lows us to conditionally select the stages of evolution of 
the system and, hence, to obtain the limits of validity 
of simple phenomenological models used for the descrip- 
tion of the complicated physical phenomenon, as well as 
to describe separate stages of development of the system 
in more detail. In particular, the stages of formation 
(t ~ 0.1 fm/c) and thermalization (r ~ 1 fm/c) are most 
often described with the use of microscopic models based 
on the processes of interaction of quarks and gluons 
3] . To describe the stage of spreading of a dense medium 
(1 < t < 7 fm/c), the relativistic hydrodynamics is most 
often in use (j-Q • The further evolution of a hadron gas 
and the process of kinetic freeze-out (r ~ 20-25 fm/c) 
are covered by kinetic models 0, Q • As parameters for 
the determination of the stages of evolution of a system, 
one can take the energy density, mean free path, rate of 
collisions of particles, etc. 

In the present work, we use the hadron reaction rate 
(number of reactions in a unit volume per unit time) in 
a given four-dimensional region of space-time as a pa- 
rameter of the spatial evolution of the interacting sys- 
tem. Such a quantitative estimate allows us to define 
the reaction zone, whose study offers the possibility of 
establishing the space-time structure of a fireball from 
the viewpoint of the interaction intensity at every point 
of space-time (for a study of the fireball structure in the 
momentum space using the number of reactions, see Refs. 
[1G3). At the same time, the regions of a fireball can be 



distinguished by the total number of collisions which oc- 
cur there. We use this quantity to determine the bound- 
aries between different zones of reactions. 

Zones of reactions. The number of reactions in the 
given space-time region can be determined with the use 
of the distribution function f(x,p). For example, in the 
approximation of two-particle reactions 2 — > 2, this func- 
tion satisfies the Boltzmann equation [11] 

Pi d»fi= [ I I W 12 ^ 34 (/ 3 h - h h) , (1) 

J2 J3 Ji 

where the right-hand side contains the collision integral. 
The quantity Wi2->34 is the transition rate which in- 
volves the reaction cross section and the conservation 
laws, J t = j ( 27r )f^. ; fi = f{x,Pi) are one-particle dis- 
tribution functions, and Ei — y/ mf + p\ is the energy 
of a particle with momentum p i . The probability of 
a collision of two particles with momenta p 1 and p 2 
corresponding to the distribution functions fi and f 2 , 
respectively, is determined at a space-time point x as 
/3J4 ^12^34 f{x,pi) f(x,p 2 ). By integrating over the 
momenta of particles, we obtain the rate or the four- 
density of reactions at the point a;: 

r(*)=/ / f fwi2^ 3i f(x, Pl )f(x,p 2 ). (2) 

Jl J2 J3 J i 

Then, the number of reactions in the given space-time 
region fi is 

JVcoii(fi)= / d A xT{x). (3) 

It is seen that the number of reactions in the given space- 
time region depends on the four-density of reactions T(x), 
which can be determined in a certain model approxima- 
tion, e.g., like that in Refs. fl2l - fl4j ]. In particular, T(x) 
can be determined with the use of transport models. 

Let us consider a large space-time region containing 
99.99 % of all two-particle reactions and decays of reso- 
nances related to the event (a particular nucleus-nucleus 
collision). That is, the separated space-time region is so 
large that the dominant part of all reactions of hadrons 
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FIG. 1: (Color online) Algorithm of 
selection of pixels for the determina- 
tion of the zone of reactions. 
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occurs in it. It can be a four-cube of reactions Cr with 
edges Li, where i — t, x, y, z. To determine the zone of 
reactions, we divide the cube into separate equal parts 
(pixels), i.e., elements of the four-space, see Fig. [T] (left 
panel). Let SI = Cl(t,r) be the four-volume of a pixel 
with coordinates of the center of this volume (t,r). To- 
tally, we have N p i x = L t L x L y L z /Q pixels. Then, for 
each four numbers (t,r), we can calculate the absolute 
number of reactions in the given pixel tt(t,r) by using, 
e.g., formula ([3]). After this, we sort the pixels from left 
to right by the following hierarchy: from a pair of pixels, 
the left pixel is that in which a larger number of reac- 
tions has occurred. The arrangement of pixels is shown 
in Fig. [1] (right panel). The total area of the whole his- 
togram (area covered by all bins) is equal to the total 
number of all hadron reactions iVtot in the four-cube of 
reactions Cr. 

It is obvious that the problem of the determination of 
a space-time region, where all collisions have occurred, 
has a statistical character. Therefore, it is reasonable to 
search for this region with a certain precision. A small 
part of unaccounted reactions, to which we can refer also 
the decays of long-lived resonances, happens on such time 
intervals and at such distances which exceed the sizes of 
the four-cube of reactions Cr. 

Let us sum the areas of bins beginning from the left 
according to the obtained hierarchy. We recall that the 
area of each bin gives the number of reactions in the cor- 
responding pixel. In such a way, we can reach the value 
of the sum which is equal to a given number (ai -/Vtot), 
where ot\ is a part of the absolute number of all reac- 
tions iVtot (see Fig. [T]). The hyperspace region which is 
occupied by the pixels contributing to this sum gives the 
reaction zone with the most intense reaction rate, and 
an ai piece of all hadronic reactions can be attributed to 
this four- volume. 

Exploiting this algorithm, we follow to sum the area 
of bins to obtain a next (a^iVtot) number of reactions. 
These bins determine a neighboring four- volume region, 
where an ct2 piece of all hadronic reactions has occurred. 
Next, we determine a four-region which contains an 0:3 
piece of all reactions, and so on. Then, a.% + ai + + 
■■■ = 1. 

In the framework of this approach (algorithm), one can 
analyze a space-time structure of the fireball for different 



species of particles and for different types of reactions, 
e.g., for elastic or inelastic reactions, decays, and so on. 
Obviously, this information gives insight into the dynam- 
ics and the structure of a particular nucleus- nucleus col- 
lision. 

Results of calculations. To carry out calculations, 
we use the transport model UrQMD v2.3 QJ, EH which 
allows one to calculate the four-density of reactions at ev- 
ery point of the space-time region and to select reactions 
of a given type and for the given species of particles. 
In the present paper, we investigate, first, all possible 
hadron reactions and, second, just inelastic hadron reac- 
tions. 

We take the number of reactions N co \\[Q,(t, r)] in a pixel 
fi(i, r) as a result of the averaging over 1000 events. In 
the calculations, we took a four-cube of reactions Cr with 
the size of edges Li = 200 fm, where i = t, x, y, z. 

In Figs. [2] and [21 we show the results of calculations for 
conditions at the BNL Alternating Gradient Synchrotron 
(AGS), Au+Au at 10.8 A GeV, and at the CERN Su- 
per Proton Synchrotron (SPS), Pb+Pb at 158A GeV, 
in the case of central collisions. In accordance with 
the proposed algorithm, we determine the four-volume 
which contains 99% of all hadronic inelastic reactions, 
2 — > 2' + m, m > [depicted as the medium-gray (red) 
area]. We name this zone as a region of hot fireball. We 
determine also a four-volume that contains 99% of all 
possible hadronic reactions which include, of course, the 
previous zone. We name the region of 99% of all hadronic 
reactions excluding the zone of the hot fireball as a cold 
fireball [dark-gray (blue) area]. In fact, one can prefer 
another precision for the determination of zones of reac- 
tions. 

By continuing to move along the hierarchical structure 
of pixels from left to right and by gathering a sum of the 
areas of their bins, i.e., the numbers of hadronic reactions 
corresponding to pixels, we determine those pixels which 
give, in total, for example, 0.9% of the total number of 
all reactions A^ot- This additional region can be named 
as a fireball halo [light-gray (cyan) area]. That is, the 
three space-time regions (hot, cold, halo) cover 99.9% of 
the total number A^ot of all hadronic reactions. 

In Figs. [2] and [31 the rate of reactions is represented 
in the coordinates z-t. That is, we deal with the cor- 
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FIG. 2: (Color online) Projection of 
the reaction zone on the z-t plane un- 
der the AGS (Au+Au at 10.8^ GeV) 
conditions. The medium-gray (red) 
region contains 99% of all inelastic 
reactions, 2 — > 2' + m, m > 0. The 
medium-gray (red) and the dark-gray 
(blue) regions together contain 99% 
of all hadronic reactions. The light- 
grey (cyan) region contains 0.9% of 
all hadron reactions only. The solid 
line bounds the region containing 
99% of all reactions. The dashed- 
dotted line bounds the region con- 
taining 90% of all reactions. 
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FIG. 3: (Color online) Same as 
Fig. [2] but for calculations under SPS 
conditions (Pb+Pb at 158 A GeV). 
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responding projection of the four-volume of the zone of 
reactions on the z-t plane. To construct this projection 
of the three-dimensional spatial pattern onto the z axis, 
we sum first all collisions along the transverse direction 
at the fixed coordinates (t, z), namely [see Eq. ||3J)], 



N co n[Sl(t,z)] = / dxdyN con [n(t,x,y,z)] 



(4) 



Thus, we put a number of reactions N co \\ in correspon- 
dence with the pixel £l(t,z) with the coordinates (i, z). 
Then we construct the hierarchy of pixels £l(t, z) accord- 
ing to the above-mentioned sorting algorithm (see Fig.Q]). 
It is seen that the interacting system has a comparatively 
long lifetime: A hot fireball decays completely only in 
the time intervals of the order of 15 fm/c for AGS and 
35 fm/c for SPS. A cold fireball lives for 30-33 fm/c for 
AGS and 90-100 fm/c for SPS. Black solid lines show the 
approximation of the boundaries of the reaction zones. 

By definition, the zones of hot and cold fireballs to- 
gether contain 99% of all reactions. Moreover, the greater 
part of reactions outside of these zones consists of decays 
of resonances. If we follow the "classical" definition of 



the sharp kinetic freeze-out hypersurface as some bound- 
ary that separates the interacting system from the space 
domain where particles do not interact, then we can re- 
gard the hypersurface which bounds the region of a cold 
fireball [the boundary between the dark-grey (blue) and 
light-grey (cyan) regions] as the sharp freeze-out hyper- 
surface. 

For both collision energies of nuclei under considera- 
tion, the upper parts of freeze-out hypersurfaces in Figs. 
[H and [3] are the space-like hyperbolas, which have the 
form of constant proper-time surfaces to within some fac- 
tor, namely, t{z) = A\> Vo 2 + z 2 , where A — 0.65, to = 
46 fm/c for AGS energies and A = 0.95, tq = 38 fm/c for 
SPS energies. From the bottom a cold fireball is bounded 
by the time-like hypersurface, which has the form of a 
straight line t{z) = in + \ z i where to is approximately 
zero, and v — 0.8 for AGS energies and v ~ 0.98 for SPS 
energies. 

The results of evaluations of the time dependence of 
the reaction rates integrated over the total spatial vol- 
ume (frequency of reactions) are depicted in Figs. 0] and 
[5] The thick solid line indicates all reaction rates in the 
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FIG. 4: (Color online) Frequency of 
the hadron reactions (reaction rates 
integrated over the total spatial vol- 
ume) for AGS conditions (Au+Au at 
10. 8 A GeV). Different curves corre- 
spond to different types of reactions. 
Vertical lines indicate the time mo- 
ments of the margins of zones of reac- 
tions for the point z — (see Fig. [2} . 




FIG. 5: (Color online) Same as 
Fig. m but for the SPS conditions 
(Pb+Pb at 158 A GeV). 



fireball, the thin solid line indicates just the elastic scat- 
tering of hadrons (2 — > 2), the dash-dotted line shows 
all inelastic reactions (2 — ^ 2' + m, where m > 0), the 
dotted line stands for fusion reactions (2 —> 1'), and the 
dashed line distinguishes decays (1 — > 2' + to, to > 0). 
Note that the compound-elastic reactions (pseudoelastic 
in Ref. [l7|) such as 2 — > p — > 2 are separated in two 
parts; the first part 2 — > p is incorporated to the reac- 
tions of fusion, and the second part p — > 2 to the reactions 
of decay. 

The main feature of the frequency of all reactions 
(thick solid lines in Figs. [4] and [5]) is its increase up to 
t « 3.9 fm/c for AGS energies and t ps 0.84 fm/c for SPS 
energies. This can be explained by increasing the number 
of nucleons as participants of the reactions, when one nu- 
cleus penetrates into another one. Indeed, the maximum 
overlap of two nuclei happens if their centers coincide. 
This time can be estimated as 



where Rq is the nucleus radius, v — pq z / \J Mfj + Pq z , 7 = 
1/y/l — v 2 , poz is the initial nucleon momentum in the 



cm. system of two nuclei, and Mn is the nucleon mass. 
In what follows, we name t c as the fireball formation time. 
For two experiments under consideration, this gives t c = 
2.9 fm/c for AGS (10. 8 A GeV) and t c = 0.77 fm/c for 
SPS (158A GeV). These values are very close to the time 
moments which correspond to the first maximum of the 
frequency of all reactions (thick solid lines in Figs. 0] and 
[5]) and local maxima of the particular types of reactions 
for SPS conditions, see Fig. [5] Slight difference of t c and 
the time point of the real maximum can be explained 
by some decrease of a nucleon velocity which is due to 
inelastic and elastic reactions (stopping) of nucleons. The 
same behavior can be seen as well in Ref. [T7j . Fig. 1, 
obtained from UrQMD calculations for SPS conditions 
(Pb-Pb, E kin = 158A GeV). 

Inelastic nucleon collisions dominate at the first stage, 
t ^ tc , of nucleus- nucleus collision. At later times, the 
elastic and fusion reactions become more significant. We 
note that for both energy conditions, the decay processes 
become the dominant ones after t ps 10 fm/c. 

After the full overlap of the nuclei, the created system 
begins to expand in space, which results in a decrease of 
the reaction rates. At the same time, the number of sec- 
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ondary particles still increases, resulting in an increase 
of the total reaction integral rate. Hence, the rate of ex- 
pansion of the system and its ratio to the creation rate of 
secondary particles will determine the result of the com- 
petition of these two tendencies. A similar competition 
was investigated in Ref. [13], where on this basis some 
conclusions about freeze-out were made. 

The number of secondary particles (mainly 7r mesons) 
is approximately (n„) sa 1.6 per nucleon for AGS condi- 
tions and (n w ) ~ 6 for SPS conditions. Such an increase 
of the number of secondary particles with the collision 
energy leads to the sufficient difference of dependences 
of the frequency of reactions on time. Namely, for AGS 
conditions after i c , the the frequency of all reactions goes 
down, which results in fireball division into two parts at 
the time moment t = ifd and the further breakup. The 
fireball division time is defined as the minimum value of 
time on the space-like hypersurface, which bounds the 
region of the cold fireball (blue area) from above, i.e., 
i fd = t(z)| (see Figs. [5] and [3]). For SPS conditions, 
one can see the second local maximum of the frequency of 
reactions (see thick solid line in Fig. [3]) at tM ~ 10.6 fm/c 
which is a consequence of a large number of reactions 
with secondary particles. This results in an increase of 
the fireball lifetime: 

t = t fd -t c . (6) 

We note that the time moment ifd depends very weakly 
on the collision energy. 

It is seen that at z — in the cm. system after the 
time moment ifd, the rates of elastic and inelastic re- 
actions vanish. That is, since this moment, the system 
behavior is determined mainly by the individual prop- 
erties of particles (basically resonances). That is why, 
in spite of the sufficient difference of collision energies of 
the experiments under consideration, the times ifd are 
approximately the same. If we compare the longitudinal 
sizes of the fireballs 2R Z at the time moment t = ifd, we 
see that they are approximately the same and equal to 
Rz — v ifd [v is defined in Eq. ((SJ), see Figs. [2] and [3]. This 
fact can explain the weak dependence of the pion inter- 
ferometric radius Rl on the beam energy, Rl oc R z [HI]. 
It can be claimed that the fireball achieves its maximum 
longitudinal size at the time moment i = ifd, when it is 
divided into two parts. 

Discussion and Conclusions. The main idea of the 
method proposed for analyzing the fireball structure is 
the determination of the hierarchy of hadronic reactions 
in accordance with their intensity (rate), see FigJTJ In 
view of the hierarchy obtained, we formulate an algo- 
rithm for the determination of the zones of reactions of 
the interacting system: we divide the fireball space-time 
volume into different regions. By definition, we recog- 
nize a specific zone of reactions ( elastic, inelastic, etc.) 
as a four-volume, where 99% of a given specific type of 
reactions have occurred. 

In the present microscopic study, we separate a fire- 
ball into the following regions, which characterize its 



evolution (see Figs. [2] and [3]): (1) a hot fireball region, 
where 99% of all inelastic hadronic reactions have oc- 
curred (medium-gray or red), (2) a cold fireball region 
(dark-gray or blue), which together with the hot fireball 
contains 99% of all hadronic reactions At t, and (3) a 
fireball halo, where 0.9% of all hadronic reactions, i.e., 
0.009At o t, have occurred (light-gray or cyan). Two last 
regions together are a space-time region containing the 
hadron-resonance gas, and the reactions in this region 
are mainly presented by decays of resonances. 

An important question that can be clarified by the 
study of the zone of reactions is how the space-time 
boundary of a fireball is related to the so-called sharp 
freeze-out hypersurface. In the literature, the sharp 
freeze-out hypersurface is defined with the help of some 
parameter P(t, r) which takes the critical value P c on 
the hypersurface. That is, the equation of the hypersur- 
face has form P(t,r) = P c . As such a parameter, one 
may choose the energy density e(i,r) 0, m , temperature 
T(t,r) dHH, density of particles n(t,r) (HI, etc. It is 
necessary to note that any definition of the sharp freeze- 
out is possible just with chosen accuracy. 

One can follow the "classical" definition: the sharp 
kinetic freeze-out hypersurface is an imaginary hyper- 
surface, outside of which there are no collisions between 
particles of the system. In other words, the sharp ki- 
netic freeze-out hypersurface is determined as a surface 
bounding the space-time region, in which almost all col- 
lisions between hadrons of the system happened (for ex- 
ample, 99%). Thus, we can identify the reaction zone 
boundary and the sharp kinetic freeze-out hypersurface. 
Then the boundary between the zone of a cold fireball 
and the fireball halo, can be interpreted as a sharp ki- 
netic freeze-out hypersurface (see Figs. [2] and [3j. In the 
coordinates (i, z), the space-time part of this hypersur- 
face is a hyperbola and has the form t(z) = A\Jto 2 + z 2 , 
where A — 0.65, To — 46 fm/c for the AGS energy 
(E kin = 10.SA GeV) and A = 0.95, r = 38 fm/c for 
the SPS energy (E kin = 158 A GeV). It is clear that the 
fireball division time is related to the parameters of the 
hyperbola in the following way ifd — Atq and we ob- 
tain ifd w 30 fm/c for AGS (E kin = 10. 8 A GeV) and 
ifd « 36 fm/c for SPS (E kin = 158A GeV). 

The lower time-like hypersurface bounding a cold fire- 
ball has the form of a straight line t(z) — i + -z, 
where to is close to zero, and v — 0.8 for AGS ener- 
gies and v = 0.98 for SPS energies. For AGS energies, 
-Ekin = 10. 8 A GeV, the time-like boundaries of the reac- 
tion zones differ significantly from one another and the 
light cone (see Fig. [2]). However, at higher SPS energies, 
for example, at E k - m — 158A GeV, the time-like hyper- 
surfaces bounding all three zones of a fireball practically 
coincide with one another and are close to the light cone 
(see Fig. |31). Thus, we can predict that this behavior will 
result in all time-like hypersurfaces merging on the ener- 
gies available at the BNL Relativistic Heavy Ion Collider 
(RHIC) and coinciding with the light cone. 

On the base of the same ideology, we can define the 
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sharp chemical freeze-out hypersurface. Namely, the hy- 
persurface separating the zones of the hot and cold fire- 
balls can be associated with the sharp chemical freeze-out 
with an accuracy of 99% (we assume that the chemical 
freeze-out occurs when the inelastic reactions are com- 
pleted [12). 

On the other hand, if one deals with a continuous 
freeze-out, for instance, by fixing the coordinates of the 
last collisions, then the set of these points (with an accu- 
racy of 99%) will be inside the boundaries we determined 
above. 

On the basis of the analysis of the time dependence of 



the frequency of reactions (see Figs.|4]and[5]), we conclude 
that there are two specific time points in the evolution of 
a hadron fireball: the fireball formation time t c defined 
as the time of the full overlap of two nuclei [Eq. ([SJ] and 
the fireball division time ifd which corresponds to the 
separation of the fireball into two individual parts. 
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